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ABSTRACT 

We report optical long-slit spectra and direct imaging (ground-based and with 
HST) of the pre-planetary nebula (pPN) M2-56 obtained at different epochs. 
The optical nebula is composed by shock-excited material distributed in two 
pairs of nested lobes with different sizes and surface brightness. The compact, 
bright inner lobes (ILs) have an angular size of ~l''5xl" each and display closed, 
bow-shaped ends. The extended, faint outer lobes (OLs), which enclose the inner 
ones, have an agular size of ~13"xl0". Within the ILs and the OLs the veloc- 
ity increases with the distance to the center, however, the ILs show expansion 
velocities larger than the OLs. Consistent with the large speeds reached by the 
ILs (of up to ~350km s~ x at the tips), we have measured the expansive proper 
motions of the knots (A^ t ~0'. / 03 yr" 1 ) by comparing HST images taken in 1998 
and 2002. Moreover, we have discovered remarkable changes with time in the 
continuum and line emission spectrum of M2-56. In 1998, we detected a burst 
of Ha emission from the nebula nucleus that is interpreted as an indication of a 
dense, fast (~350-500km s -1 ) bipolar wind from the nebula's core (referred to 
as "Fl-wind"). Such a wind has been recently ejected (after 1989) probably as 
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short-duration mass- loss event. Our data also reveal an optically thick compact 
cocoon (or shell-like structure) and a H II region around the central star that 
result from further post-AGB mass-loss after the Fl-wind. Recent brightening of 
the scattered stellar continuum as well as an increase of scattered Ha emission 
along the lobes is reported, both results pointing to a decrease of the optical 
depth of the circumstellar material enshrouding the star. The data presented 
here unveil the complex post-AGB mass- loss history of this object, whose rapid 
evolution is driven by multiple episodes of mass outflow, not regularly spaced in 
time, leading to: (i) acceleration of the molecular envelope that surrounds the 
optical nebula (kinematical age t k ~1400yr - Castro- Carrizo et al., 2002), (ii) the 
OLs (t k ~350-400yr), (in) the ILs (t k ~40yr), (iv) the Fl-wind (t k <10yr), and 
(v) the nuclear cocoon and H II region (t k <2yr?). The sucessive multiple post- 
AGB winds in M 2-56 are characterized by ejection speeds increasing with time. 
In contrast, the mass-loss rate and linear momentum show a time decreasing 
trend. 

Subject headings: stars: AGB and post-AGB, stars: mass loss, circumstellar 
matter, ISM: jets and outflows, planetary nebulae: general 



Introduction 



Planetary nebulae (PNe) evolve from Asymptotic Giant Branch (AGB) stars after a brief 
(«10 3 yr), intermediate stage known as the pre-Planetary Nebula (pPN) or post-AGB phase. 
During the AGB to PN transition dramatic transformations occur: the spherical, slowly ex- 
panding (Vexp ~ 15 km/s) circumstellar envelope (CSE) expelled during the AGB becomes a 
nebula with, usually, clear departures from sphericity and fast (~ 100 km/s) outflows directed 
along one or more axes. The varied PN and pPN morphologies include not only axisym- 
metric (elliptical and bipolar) shells but also multipolar (multiaxial) structures, multiple 
co-axial shells, hig hly collimated jet- l ike ejections sometimes arranged in a point-s ymmetric 
structure, etc (e.g. ISahai et al.lll998l . hQQ'A lUeta et al.ll2007l : Isiodmiak et~aDl2008h . Most of 
the observed PN and pPN morphologies are difficult to explain without invoking; the action 



of processing, perhaps epis odic, collimated outflows (e.g. iMiranda et al.lll999t I Guerrero et al. 
200 ll : IMiranda et al.ll2006l ). In fact, collimated outflows are conside red as the main shapin g 
agent of PNe through their interaction with the spherical AGB CSE ( ISahai fc Traugerlll998l ). 
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The existence of very young, high velocity bipolar jets in some P PNe and evolved AGB 
stars is evidenced, for e xample, by water maser observations (e.g. Ilmai et al.l l2002r llmai 



20071 ; ISuarez et al.ll2009l . and references therein). In some pPNs/PNs, fast, post-AGB winds 
are also revealed by P-C ygni profiles in the Ha (and other recombination lines) emissio n 
close to the central star ( iSanchez Contreras fc Sahaill200ll ; ISanchez Contreras et al.ll2008l ). 
In spite of the growing evidence of jet-like ejections in PPNs, their origin, nature (episodic 
or continuous?), typical mass- loss rates, life-times, etc, are still very poorly known. Detailed 
studies of PPNs are crucial to obtain information about the properties and dynamics of 
post-AGB winds and their evolution. Optical spectroscopic observations of pPNs and PNs 
are particularly useful for probing the post-AGB winds and their interaction with the CSE 
formed in the previous AGB phase. 

This work focuses on the pPN M2-56 (= IRAS 23541+7031= PK 118+08), which so 
far remains relatively poorly characterized, especially in the visible. A ground-based optical 
image obtained by Goodrich (1991, hereafter G91) showed a ~10"x2" bipolar nebula with its 
long axis oriented in the East- West direction (P. A. =90°). The limited quality of this image 
(PSF~l'/6 and low S/N) impeded a correct identification of the different nebular components 
in M2-56 revealed by data presented in this work. High-angula r resolution opt i cal im ages 
obtained with the HST hav e been used as complementary data by lSiodmiak et al.l (120081 ) and 
Castro-Carrizo et al.l ( 120021 ) although a detailed study based on these images has not been 
done so far. The HST images show that the bright nebulosity referred to as the "West lobe" 
by G91 is resolved into two compact lobes (~l'.'5x 1") emerging from the center. The emission 
line optical spectrum of M 2- 56, with strong forbidden lines of low-excitation ions, in dicates 
major excitation by shocks ( )Goodrichlll99ll ; iTrammell et al.lll993l ; iRiera et al.l 120061 ) . The 
lack of substantial polarization in the observed line spectrum implies that the latter is locally 
produced in the shock-excited lobes of M 2-56 and that there is not a significant contribution 
by line emission from the central regions reflected by the nebular dust (" scattered spectrum" , 
Trammell et al.lll993l ). (In this work we show that this situation has recently changed). The 
optical spectrum of M2-56 also shows a weak red continuum that is the l ight from a B- 
type central star attenuated and/or partially reflected by circumstellar dust (I Cohen &: Kuhi 
19771 1 . 



A detailed study of the dynamics and mo rphology of the molecular en velope of M 2-56 
with high- angular resolution was carried out by lCastro-Carrizo et al.l ( 120021 ). Interferometric 
CO mapping shows an incomplete ~30"xl5" hourglass-like envelope surrounding the optical 
lobes and similarly oriented to the latter in the plane of the sky. The inclination of the 
symmetry axis of the CO envelope with respect to the plane of the sky is estimated to 
be z=17°±2°, with the East side of the nebula tilted towards us. The densest molecular 
material lies at the center, where the intersection of the two bipolar lobes shapes an equatorial 
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ring/torus expanding at 7-8 km s -1 . Two compact CO clumps are found at the tips of 
each lobe expanding axially with velocities of ~100-200km s^ 1 . The total mass and linear 
momentum of the molecular envelope of M2-56 is M mo i=0.05 M & and P mo i=2M Q km s -1 , 
adopting a distance to the source of d—2.1 kpc and a standard value of the CO-to-H 2 relative 
abundance of X(C0)=2x 10~ 4 . The morphology and dynamics of the bipolar nebula is 
interpreted as the result of the interaction between fast, collimated post-AGB ejections and 
the slow CSE formed during the AGB phase. The Hubble-type expansion of the molecular 
outflow suggests that such an interaction took place approximately 1400 yr ago in a relatively 
short time, <300yr. 

In this paper we will ado pt a distance to M2- 5 6 of rf=2.1kpc based on the detailed 



discussion and final estimate by lCastro-Carrizo et al.l (120021 ). As estimated by these authors 



for such a distance, the luminosity of M2-56 obtained by integrating its spectral energy 
distribution (SED) from the optical to the mm-wavelength range is ~55OOL . 

In this work, we report multi-epoch optical long-slit spectra and images of M 2-56 that 
have allowed us to describe the structure, kinematics, physical conditions, and mass-loss 
history of this object with unprecedented detail. Our results are discussed in terms of the 
formation and evolution of M2-56, which is currently in rapid transition from the AGB to 
the PN stage and presents clear signs of on-going, variable wind activity at the nucleus. A 
summary of our findings is given in Section [HJ 



2. Observations and data reduction 
2.1. Optical Imaging 

2.1.1. Ground-based images 

We have obtained broad- and narrow-band images of M2-56 using the Wide Filed Cam- 
era (WFC) at the prime focus of the 2.5m INT of the Roque de los Muchachos Observatory 
(La Palma, Spain). Observations were performed on service time mode on June 13 , 2009. 
The WFC consists of 4 thinned EEV 2kx4k CCDs. The CCDs have a pixel size of 13.5 /im 
corresponding to 0'.'33 per pixel. Our target was observed with the chip CCD4, which has a 
field of view of ~12'x23'. Two filters were used: Ha (A C =6568A, FWHM=95A) and Harris 
R (A C =6380A, FWHM=1520A). The integration times were about 530 and 2400 s for the R 
and Ha images, respectively. The compact, nuclear regions of M2-56 are saturated in the 
Ha image. The weather was good during the observations and the seeing was about l'.'S-l'.^. 

Images have been debiased and flat field corrected following the standard procedure 
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using irafB Individual frames have been registered and combined for each filter to produce 
final images shown in Fig. [I] 

2.1.2. HST images 

We have used two (Ha and continuum) high angular resolution images of M 2-56 from 
the HST archive. The Ha image was obtained with the WFPC2 PC camera using the narrow 
band filter F656N (A C =6564A, FWHM=21.5A) on 1998-09-04 (GO 6761; PI: S. Trammell). 
This camera has a field of view (FoV) of 40"x40" and a plate scale of 0'.'05 per pixel. Total 
exposure time was 5000 s. The continuum (V band) image was obtained with the ACS/HRC 
camera using the broad-band filter F606W (A C =5907A, FWHM=2342A) on 2002-12-23 (GO 
9463; PI: R. Sahai). The ACS camera has a field of view of 26"x29" and a plate scale of 
0'/025 per pixel. Total exposure time was 800 s. The images (Figs. [I] and |2]) were pipeline 
reduced. 



2.2. Optical Spectroscopy 

We have obtained optical long-slit spectra of M 2-56 using the Intermediate Dispersion 
Spectrograph (IDS) of the 2.5m Isaac Newton Telescope (INT) and the Andalucia Faint 
Object Spectrograph and Camera (ALFOSC) mounted on the 2.6 m Nordic Optical Telescope 
(NOT) of the Roque de los Muchachos Observatory (La Palma, Spain). Observations were 
performed in three different campaings in 1998, 2000, and 2009. The log of our multi-epoch 
spectroscopic observations is given in Table[TJ 

2.2.1. INT+IDS 

For the first set of spectroscopic observations with IDS (run#l), the detector used was 
a TEK 24 /im CCD with 1124x1124 pixels mounted on the 500mm camera. Three slits 
positions were observed: one at position angle PA=90° (the symmetry axis of the nebula) 
passing through the nebula center (referred to as PA90) and two oriented along PA=0°, one 
passing through the bright nebula center (PA0) and the other approximately through the 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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center of the East lobe (at offset Aa~+6" from the center, PAOE) - see Fig.[3l leftmost 
panel. The slit was 1" wide and long enough to cover the whole nebula. CuNe lamps were 
used for wavelength calibration. The achieved spectral resolution (FWHM of the calibration 
lamp lines) was ~44km s" 1 around Ha. 

For the second set of spectroscopic observations with IDS (run#2) the detector was 
the EEV10 CCD, with squared pixels of 13.5/im lateral size. Only a clear and unvignetted 
region of 700x2600 pixels was used (the 2600 pixels were along the spectral axis). The CCD 
was mounted on the 500mm camera. A total of four slits positions were observed with the 
R1200Y grating (Fig.[3l leftmost panel): we repeated positions PA90 and PAOE, already 
observed in run^l, and also added two more, one oriented along PA=0° passing through 
the West lobe (at offset Aa~- 6" from the center, PA0W) and another one passing through 
the nebula center but oriented along PA=119° (PA 11 9). For the slit position PA90 we also 
observed with the R900V grating. We used CuNe and CuAr lamps to perform the wavelength 
calibration. The spectral resolution achieved was ~50km s -1 , at Ha, and ~83km s" 1 , at 
H/3. 

In both runs, #1 and #2, flux calibration was done using sensitivity functions derived 
from spectrophotometric standards, namely, HR1544, HR3454, and HD217086, and taking 
into account the atmospheric extinction curve at La Palma. All observations were performed 
under photometric conditions and with a typical seeing of ~1" except for those carried 
out the first night of our run ^2 (2000-11-11), during which an average seeing of ~2" 
was measured. A secondary flux calibration of the mentioned dataset observed under non- 
photometric conditions (i.e. within the ~6200-6800A range covered by grating R1200Y) was 
performed by matching the fluxes of the [S II]AA6716,673lA doublet lines measured in 1998 
and 2000, that is, assuming that the line fluxes have not changed with time. This is a 
reasonable assumption given the stability of the profile of these lines in all our campaigns 
(including run #3 in 2009). We derive a value for the scale factor of ~4 (i^ooo x 4~^i998)- 
Once applied to all line fluxes in the range 6300-6740A, this factor also corrects for the 
different seeing in our non-photometric data (~2"), which affects the fraction of the nebular 
emission that enters the slit. We double-checked the obtained flux scale factor by comparing 
the fluxes integrated over the nebular core (the inner ~2") of all lines (not only [S II]) in 
the spectra taken in 2000 along PA90 (under non-photometric conditions) and along PA119 
(under photometric conditions), since both slits pass through almost the same central regions. 

The data were reduced using standard procedures for long-slit spectroscopy within the 
IRAF package. The long-slit spectra presented in Figs.|3] and H] have been smoothed with a 
flat-topped rectangular kernel of dimension 3x3 pixels to increase the S/N. 
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2.2.2. NOT+ALFOSC 

Observations were done in service time mode using the CCD #8 detector, with 2048 x 2048 
13.5 fim sized pixels, and the VPH Grism #17. A O'.'9-wide (5(3- long) slit oriented at PA=90° 
and passing through the nebula center (slit PA90) was employed. The effective spatial res- 
olution (seeing) was 0'.'9. We observed and used Ne lamps for wavelength calibration. The 
achieved spectral resolution ~61 km s" 1 around Ha. The spectrophotometric standard star 
BD +2842 11 was observed and used for primary flux calibration. A secondary flux scale 
factor of ~7 has been obtained, following the procedure described in the previous subsection 
12. 2. 1| and applied to the data. Final spectra, obtained after standard data reduction with 
IRAF, are shown in Figs. 013 

3. Results 
3.1. Imaging 

In Fig. [T] we present images of M2-56 that unveil the different components of the optical 
nebula and their morphology. We find two pairs of nested lobes similarly oriented in the 
plane of the sky (along P.A.~90°) but with different sizes and a large contrast in surface 
brightness. The bright inner lobes (ILs), which can only be identified and spatially resolved 
in the HST images, have an angular size of ~l'.'5xl" and display closed, bow-like ends, 
which very likely represent shock fronts. A 'three finger'-like structure at the base of the 
east IL (within 0'.'5 from the core) can be seen in the F606W image. The faint, more extended 
(~13"xl0"), outer lobes (OLs), which are probed by our deep ground-based images but are 
only barely detected in the HST images, enclose the inner ones. Both pairs of lobes present 
an overall axial symmetry, however, they display a remarkable clumpy structure, including 
strings of compact knots and curved features that do not always have an axial counterpart. 

The surface brightness of the east lobes is larger than that of the west onefl This is 
most likely a result of a larger extinction of the west side of the nebula, which is receding 
from us, given the hourglass distribution of the molecular, and presumably dust, envelope 
around the optical nebula (§[T]and §EJ). The larger brigthness contrast at the central parts of 
the nebula (e.g., note that the three-finger structure of the elL does not have a counterpart 
in the west) is consistent with a larger extinction along the densest equatorial regions. 

In spite of their different filter bandwidth, both the broad- and narrow-band HST 



2 Hereafter, we will refer to the east (west) ILs and OLs as elLs and eOLs (wILs and wOLS), respectively. 
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images of M2-56 trace approximately the same gas component, as deduced from our long- 
slit spectra (§ 13. 2p . First, the contribution of the continuum to the emission within the 
F606W and F656N filters was negligible at the time the HST images were acquired, except 
maybe in the innermost regions of the ILs (within ±0'.'25) where a very weak continuum level 
is detected in our spectra. Second, although the broad-band filter F606W includes several 
lines other than Ha, e.g., [N I]A5200A, Na IA5893A, [O I]AA6300,6363A, [N II]AA6548,6583A 
and [S II]AA6717,673lA (see also G91), the majority of these low-excitation transitions are 
expected to arise in the same Ha-emitting regions probed by the F656N filter. The only 
exception could be the [N II] and [S II] doublets, which may proceed from regions with 
a slightly higher excitation, however, these lines have a contribution to the total line flux 
within the F606W filter of less than 30% (see Tableland G91). 

We have noticed some differences between the two HST images, F656N and F606W, 
which were observed in 1998 and 2002, respectively. In particular, the string of knots/clumps 
at the tips of the ILs in 2002 is located s^O'/lO-O'/lG ahead the same structure as it was 
observed in 1998 (Tig. [2]). This is true for both the east and west ILs and we believe it is due 
to proper motionqj (A0 t ~O / /O3yr _1 ) resulting from the fast expansion of the ILs (see § 13.2.11 
and EE} . 

The base of the elL also has a quite different appearance in the F606W and F656N 
images. In particular, the three finger-like structure visible in the F606W image is absent 
in the F656N image. In the latter image, the same region, i.e. the base of the elL, shows 
an asymmetrical cone-like morphology but has a smoother brightness distribution slightly 
elongated to the north-east. These changes cannot be unambiguously attributed to real 
variations with time of the nebular structure since we cannot rule out some contribution 
from continuum emission (due to dust scattering in these innermost regions) that may be 
different in the F656N and F606W filters. 

Finally, we have compared our two-epoch Ha images observed with the HST in 1998 
and with the INT in 2009. There is an emission clump located at offset Aa=5", A5=8" that 
is visible in both data sets (clump CI; Fig.[TJ). Surprisingly, there is another pair of even 
brighter clumps in the ground-based image, labeled C2 and C3, which are not detected in 
the HST images. Other fainter clumps or knots within the OLs (e.g. C3') visible in our 2009 
images may have remained below the noise level in the HST images. 



3 Here the term "proper motion" is not used to describe the apparent motion of the nebula as a whole 
with respect to more distant stars but to refer to the tangential motion of material in the expanding nebula 
(e.g. the knots). 
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3.2. Spectroscopy 

In this section, we show the results from our multi-epoch spectroscopic observations. 
In Table[2] we report the fluxes (F\) of the different lines detected. These fluxes have been 
measured from the PA90 spectrum integrating along the l"-wide slit and over the full width 
of the spectral profile, after fitting and subtracting the continuum level. In the table, we also 
provide dereddened fluxes for a mean logarithmic extinction coefficient of c(H/3)=1.5±0.1 
(see §EJ. 

We have discovered significant changes with time in the line and continuum emission 
spectrum of M 2-56. Amongst the different lines detected, Ha is the transition displaying the 
most extreme profile variations. We will first report our results obtained from the spectra 
observed in 1998 and 2000, for which the differences found only affect the line emission from 
the nuclear regions. Next, we will describe the long-slit spectrum along PA90 obtained in 
2009, which shows additional changes with respect to the data obtained in our two previous 
runs. 



3.2.1. 1998 and 2000 observations 

In Figs. [3] and HI we show the optical spectrum of M2-56, which is composed of recombi- 
nation and forbidden line emission plus a weak, red continuum. In 1998, the continuum emis- 
sion arises from a central compact region that is unresolved in our data (with a PSF~0'/9). 
The continuum intensity near Ha (around 6510A) in 1998 integrated over the l"xl" central 
region of M2-56 is ~1.5(±0.4)xl0 _17 ergs~ 1 cm" 2 A" 1 . In 2000, we derive an upper limit for 
the continuum intensity within the same central region of <2.4xl0 _17 ergs -1 cm" 2 A" 1 (i.e. 
<2cr). Therefore, the non detection of the continuum in 2000 is consistent with its intensity 
remaining constant from 1998 to 2000. (As shown in § I3.2.2[ the continuum is brighter in 
2009). 

The long-slit Ha spectrum is characterized by a wide profile with a quite complex 
spectral and spatial distribution (Fig.[3]). Along the nebula axis (slit PA90), the Ha emission 
can be traced from the brightest central regions up to the fainter tips of the OLs at ±12". The 
Ha emission from the east and west lobes (approaching to and receding from us, respectively, 
§[!]) is, on average, blue- and red-shifted, consistent with an expansive nebular kinematics. 
The spatial distribution of the emission from the bright, compact ILs is only partially resolved 
in our long-slit spectra and can be traced from the center up to ±2". In these innermost 
regions, the much fainter OLs are not expected to contribute significantly to the Ha emission 
given the large brightness contrast with the ILs evidenced by the HST images. Along the ILs, 
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the line-of-sight or radial expansion velocity increases with the distance reaching a maximum 
value of K.=Vlsr— V sys ~110km s" 1 at the tips. In these outermost regions of the ILs, the 
Ha profile is rather broad, with a full width at zero intensity (FWZI) of ~360km s _1 . The 
Ha long-slit profile along and across the OLs (e.g., along PAOE) is roughly consistent with 
the lobes being expanding bubble-like structures. We measure a value for the mean radial 
expansion velocity of the OLs of K~50km s _1 , which is smaller than that derived for the 
ILs. In the outermost parts of the OLs, at axial offsets ±10"-12" (PA90), our spectra indicate 
a slightly smaller radial expansion velocity of VJ.=30-35km s -1 . 

We have found remarkable differences between the Ha emission from the nucleus of 
M2-56 as observed in 1998 and in 2000 (Fig. [3]). In 1998: a) the Ha profile shows emission 
wings that are much broader (FWZI~1000km s" 1 ) than in 2000 (FWZI~450-500km s" 1 ) 
- this difference is not explained by the lower S/N achieved in 2000 (note that the fainter 
emission from the OLs is well probed in our 2000 spectrum); and b) the Ha emission peak is 
blue-shifted by ~110km s _1 from the systemic velocity, unlike the spectrum in 2000, which 
peaks at V^ ys . We refer to the intense, blue-shifted emission feature and its broad wings 
observed in 1998 as feature Fl. We have checked that feature Fl was not present either in 
the long-slit spectrum along the nebula axis acquired in 1989 by G91. The broad wings of 
feature Fl are not symmetric, the blue wing seems to be depressed relative to the red one, 
which may be interpreted as a blue-shifted (P-Cygni like) absorption in the profile. The 
blue-shift of the wing absorption is found at V^~— 430km s" 1 . The presence of the spectral 
component Fl is clearly appreciated also in the Ha spectrum taken with the slit oriented 
along the nebula equator PA0. The angular size of the region that produces feature Fl 
is very small, ~0'.'5 along the main symmetry axis of the nebula and unresolved along the 
equator (<0'.'6). 

In 1998, feature Fl was not only visible in the Ha spectrum but also in the [N Il]AA6548,6584A 
lines. In contrast, feature Fl is not present in the [S II]AA6716,6732A transitions: note that 
the two lines of the doublet have comparable profiles in 1998 and 2000 (see Fig. [7]). The 
profiles of the [S II] lines are indeed very much alike to those of Ha and the rest of the 
transitions as observed in 2000 (and in 2009, see next subsection § I3.2.2[) . 

3.2.2. 2009 observations 

Figs. [5] and [6] display the long-slit spectrum of M2-56 obtained in 2009. Multi-epoch 
spectra extracted over the central l"xl" region are compared in Fig.0 

In 2009, a relatively faint, red continuum is observed not only at the center (as in 1998) 
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but also along the lobes of M2-56. As for most pPNs, the continuum emission is dominated 
by the stellar photospheric continuum scattered off by the nebular dust (we have checked that 
the contribution by nebular continuum is negligible given the intensity of the Ha emission). 
The spatial distribution of the scattered continuum peaks at the nebula center and has two 
relative maxima at axial offsets ±6'.'5. The latter regions correspond to condensations C3 and 
C3' in the ground-based images (Fig. [I]). The central region where the continuum is observed 
has a deconvolved angular size along the axis of ~l'/5, i.e. significantly larger than in 1998 
(<0'.'5). The continuum flux near Ha (around 6510A) within the l"xl" central region in 
2009 is ~3.2xl0 -16 ergs -1 cm -2 A -1 , that is, larger than that measured in 1998 by a factor 
~20. 

The shape of the Ha long-slit profile shows marked differences with respect to that 
observed earlier, not only in the central regions but also along the extended OLs. The Ha 
emission from both the approaching and receding (west and east, respectively) lobes is on 
average red-shifted, in contrast to what we observed in 1998 and 2000. The red-shift of the 
Ha emission from the aproaching eOL is smaller than that from the receding wOL. At the 
nebula center, the Ha emission has a broad asymmetric profile that peaks at VLSR=65km s -1 , 
i.e. red-shifted (by ~90km s -1 ) with respect to V^ ys . The FWZI of the Ha nuclear profile 
cannot be accurately determined because of blending with the adjacent [N II] doublet; we 
estimate that it is larger than 1300km s _1 and probably reaches up to 2500km s _1 . 

The observed red-shift of the Ha line for both the approaching and receding lobes in- 
dicates that in 2009 a significant fraction of the Ha emission in the lobes is scattered, i.e. 
is not locally produced in the lobes but rather arises at the nebular core and is reflected 
by the nebular dust. In fact, the similar spatial distribution of the Ha emission and the 
reflected continuum beyond ±6" indicates that in these regions the Ha emission is mainly 
scattered. In the inner regions, however, some contribution of unscattered Ha emission 
(produced locally in the lobes) can still be appreciated. Composite Ha line profiles con- 
sisting of two co mponents (scattered a nd unsc attered) a re also fou nd in other pPNs, for 



example, Ml-92 flTrammell et al.lll993f) . M2-9 flSolJl2000h . CRL618 flTrammell et all 11993 
Sanchez Contreras et al.ll2002l ). etc. 

In contrast to Ha, the long-slit profiles of the forbidden transitions [O I]A6363A, the 
[N II]AA6545,6584A and [S II]AA6716,673lA in 2009 are totally consistent with those observed 
in our previous run^|: roughly point-symmetric with respect to the spatial origin and V sys 
and with the emission from the east and west lobes blue- and red-shifted, respectively. This 
is consistent with the forbidden lines being predominantly produced locally in the shocked 



4 Except for the presence of feature Fl in the [N II] doublet lines in 1998, § 13.2. II 
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nebular material. 

We have discovered a number of emission lines in this work, including several Fe II, 
Fe I, and Si II permitted transitions (Fig. and Table[2]). We have also detected a broad 
(FWHM~6A) absorption feature centered near 6673A at the nebula center but slightly 
blue-shifted, near 6670A, along the east and west lobes. The relative blue-shift of the band 
between the center and the lobes indicates that it does not have an interstellar or telluric 
origin but rather is associated to M2-56. We have not found a satisfactory identification 
for this broad absorption, which may be a solid state or ice feature given its large width; a 
blend of atomic lines cannot be ruled out either. 

The emission distribution of the recombination lines discovered is quite different from 
that of the forbidden transitions. As for Ha, the former are red-shifted with respect to l^ ys 
along the approaching eOL as well as at the nebula center (these transitions are not detected 
in the fainter west lobes). The red-shift measured in the eOL is smaller than that in the 
center. We conclude that the observed recombination lines most likely arise at the nebula 
nucleus and their emission is scattered by circumstellar dust. 

Another important change detected in 2009 is the increase of the relative intensity of 
Ha with respect to the various forbidden transitions observed. For example, in 2009, the 
Ha-to-[N II]A6584A intensity ratio reaches a value of ~15, much larger than in earlier epochs, 
~3.2-3.4. The large Ha-to-[N II]A6584A intensity ratio is consistent with the presence of a 
compact H II region at the heart of M 2-56 responsible for the extra Ha flux in 2009: note 
the brightnening of the Ha line by a factor RilO relative to the flux measured in 2000. 

4. Nebular components and their spatio-kinematic structure 

In this section we describe the main nebular components in M 2-56 revealed by our data 
and analyze their spatio-kinematical structure. We also estimate the inclination of the nebula 
axis, which is a critical parameter to correct for projection effects both in the geometry and 
expansion velocities observed. 

4.1. The inclination of the nebula 

The inclination of the nebula axis is derived in two different ways. First, from the proper 
motions observed in the HST images (Fig.|2]and § !3.ip . For a distance to M2-56 of <i=2.1 kpc 
(§[1]), the observed proper motions (A^t^O'.'OS yr _1 ) imply a tangential expansion velocity (i.e. 
projected on the plane of the sky) of Vt~325km s -1 . Since the line-of-sight or radial velocity 
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at the tips of the lobes measured from our long-slit spectra is V^llOkm s _1 (§ 13. 2p . we 
obtain an inclination angle to the plane of the sky of z=tan _1 (y r / Vt)=19°. This value is in 
very good agreement wi th the independent es t imate of i=17° for the symmetry axis of the 
bipolar CO envelope by ICastro-Carrizo et al.l ( 120021 ). We derive an absolute (deprojected) 
expansion velocity at the tips of the ILs of ~350km s _1 . This value is similar to the FWZI 
of Ha emission at the ends of the ILs; the latter is known to be a reliable indicator of the 
shock velocity for bow-type shocks, indepen dent of orientation an gle, preshock density, bow 
shock shape, and preshock ionization stage ( jHartigan et al.lll987l ). 



An independent estimate of % can be obtained from the radial velocities measured at 
clumps C3 and C3' in the east and west OLs, respectively, where the Ha emission is mainly 
scattered (§ !3.2.2p . Assuming that regions C3, C3', and the nebula center are aligned along 
a given axis, which is inclined to the plane of the sky by i, and that the clumps move away 
from the central star at the same expansion velocity, V exp , then we can readily find V^ xp and 
% from the following expressions: 



cxp 



sum 



V T (C3) + V T {C3') 

" 2 
[V T (C3>) - V r (C3)] 

[V X {C3') + V T {C3)\ 



(1) 
(2) 



(see, for example, ISolJ 120001 ) . Taking into account the values of VJ.(C3)=31.3 km s _1 and 
\4(C3')=60.9km s _1 measured on the Ha spectrum along PA90 (Fig. [6]), we obtain V^, xp =46 km s" 
and 2=18.8°. The value of i found is in excellent agreement with that computed for the ILs 
from the proper motions and that of the bipolar molecular envelope. 



4.2. Spatio-kinematic model of OLs and ILs 



In order to investigate the structure and velocity field of the optical lobes of M2-56 
we have developed a simple spatio-kinematical model that reproduces the 2D profile of our 
long-slit spectra. Although the model has been mainly applied to the PA90 spectrum, we 
have checked that it also qualitatively explains the spectra at the other slit orientations. 
Our model calculates the position-velocity (p-v) diagram assuming a certain morphology for 
the envelope (compatible with our direct images) and a given velocity field, and also taking 
into account the inclination of the nebula to the plane of the sky (i=18°). We note that our 
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code is not an emission model, that is, does not predict the Ha brightness distribution along 
PA90. 

In Fig.|H]we present the geometry adopted and synthetic p-v diagrams predicted by our 
best model for the Ha line. The nebular component responsible for the spectral feature 
Fl is not included in our model because it is not spatially resolved in our spectra. We are 
assuming in our model that the lobes are hollow and show an overall axial symmetry, as 
suggested by the optical appearance of M 2-56. For the OLs and the ILs, the model fit yields 
a linear size of about (4xl0 17 cm)x(3xl0 17 cm) and (6xl0 16 cm) x(4xl0 16 cm), respectively, 
measured along the symmetry axis and the perpendicular direction where the lobes reach 
their maximum extension. 

For the adopted geometry, the velocity field that best reproduces the long-slit spectrum 
of the OLs consists of three main components: 1) radial expansion with the modulus of the 
velocity at each point increasing linearly with the distance to the nebula center; 2) axial 
expansion with the velocity modulus increasing linearly with the distance to the equatorial 
plane; and 3) shear-flow kinematics, with the velocity vectors tangential to the curved surface 
adopted for the lobes and with a constant velocity modulus. The components that dominate 
the velocity field in M2-56 vary from the base of the OLs, where components 1) and 3) 
prevail, to the external regions, where the velocity field is essentially axial. The velocity field 
described above, with the flow speed increasing with the distance to the center, works fine for 
most parts in the OLs, however, is not able to reproduce the relatively low radial velocities 
measured at the tips, ~110km s -1 around ±10-12", which implies an abrupt slow down of 
the flow. In these outer regions, represented in our model as bow-shaped caps detached from 
the main body of the OLs, the expansion seems to be dominated by a constant velocity axial 
field. 

Since the spatial structure of the ILs is only partially resolved in our 2D spectra, we have 
adopted a simple axial velocity field that is able to reproduce the overall velocity gradient 
observed within the inner ±2". We deduce a rather steep variation of the deprojected 
expansion speed, increasing from tens of km s _1 near the center to ~350km s _1 at the lobe 
tips. 

4.3. The wind responsible for the transient Ha feature Fl 

We believe that the Ha emission burst observed in 1998, referred to as feature Fl, 
is the consecuence of a fast, dense, and short-lived post-AGB wind recently ejected (§E]). 
We have measured a total angular size along the nebula axis of the region emitting feature 
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Fl of ~0'.'5 (1.57xl0 16 cm). The spectral profile of Fl, notably asymmetrical, peaks at 
V r ~110kms _1 , i.e. bluewards of V^, ys . This is consistent with a compact bipolar flow 
oriented similarly to the ILs and OLs, the red-shifted emission from the receding (west) side 
of the Fl-wind being significantly extinguished by the optically thick equatorial regions of 
M 2-56. A larger extinction towards the inner regions of the wIL (compared to the base of 
the elL) is, in fact, inferred from the brightness contrast between these two regions observed 
in the HST images and is also expected from the equatorially enhanced density distribution 
of the molecular envelope that surrounds the optical nebula of M2-56 (§[T]and §E]). (A one- 
sided collimated Fl-wind is not impossible but improbable given the bipolar morphology, 
approximately symmetric with respect to to the equator, of the optical and CO envelopes 
of M2-56). Adopting a bipolar flow with an inclination to the plane of the sky of i=18°, 
we obtain a linear extension for this component of /fi~1-65x10 16 cm. The blue-shift of the 
intensity peak of Fl indicates that the bulk of the material in the Fl-wind expands at a 
mean velocity of V ex p~110/sin(18 o )~360km s , very similar to the expansion speed at the 
tips of the ILs. The broad wings of feature Fl and the blue-shift (P-cygni like) absorption 
of the nuclear Ha profile indicate smaller amounts of gas expanding at larger velocities (up 
to ~500km s _1 ) in the Fl-wind. 



4.4. The central H II region 

The presence of a compact H II region around the central star of M2-56 is inferred 
from our data (§ I3.2.2[) . In 2009, a fraction of the Ha emission from such a region escapes 
the stellar vecinity and reaches the ILs and OLs where it is scattered off by dust. The 
kinematics of the H II region cannot be straightforwardly determined from the Ha line 
profile since the emission from this region is seen only after being reflected by the nebular 
dust, the resulting line profile being significantly affected by the distribution and kinematics 
of the dust. The broad profile of the Ha emission from the nebula nucleus (FWZI>1900- 
2500km s~ x ) may indicate rapidly outflowing material in the stellar neighborhood, with very 
large expansion velocities of up to ^1200 km s^ 1 , however, we cannot rule out other line 
broadening mechanisms (e.g., Ramman scatte ring) being totally or partially responsible for 



the large FWZI of the Ha transition (see e.g. ISanchez Contreras et al.ll2008l . and references 
therein) . 

The Ha scattered profile observed at C3' and C3 represents the profile emitted by the 
H II region as seen from the dust clumps located at those positions, i.e. from a "pole-on" 
view. As we can see in Fig.[6l the Ha profiles at C3 and C3' are very similar (except for an 
overall relative red-shift of 30-40 km s _1 ): they have a narrow intense core plus broad wings 
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and show a notable P-cygni like profile, with the blue side depressed. The Ha emission 
towards the center exposes the profile emitted by the H II region as seen along the line-of- 
sight, i.e. from a direction nearly perpendicular to the lobes. Its overall red-shift indicates 
a significant contribution of scattered Ha emission also towards the nebula center. In this 
direction, the Ha profile is broader than that observed towards clumps C3 and C3' and does 
not exhibit prominent blue-shifted absortion. This suggests that the P-cygni like absorption 
observed at C3 and C3' is produced by neutral or partially neutral material outside the H II 
region and very likely distributed in a bipolar flow. 



5. Extinction 



We have used the Ha to H/3 flux ratio measured in 2000 (Table[2]) to estimate the ex- 
tinction in M2-56. Since the Ha and H/3 fluxes have been obtained integrating spatially 
the spectrum along slit PA90 over the whole nebula, the value derived represents an average 
of the extinction along the lobes, which is expected to decrease from the center (nebular 
equator) to the tips of the lobes. We have measured a Balmer decrement of Ha/H/3~12±1, 
for which we obtain a logarithmic extinction coefficient c(H/3)=1.5±0.1 adopting in intrinsic 
ratio Ha/H/3~3.8. The latter value of the intrinsic Ha/H/5 ratio is expected for the shock- 
excitation conditions in M 2-56 derived by several authors comparing obser vational line rel- 



ative intensities with t he predictions of theoretical shock excitation models (lGoodrichlll991 



Trammell et al.l Il993l ; iRiera et al.l 12006k and references therein). (For an intrinsic Balmer 
decrement closer to the recombination value, Ha/H/3~3, we obtain c(H/3)=l. 7-2.0.) From 
the value of c(H/3), we deduce an average optical depth near Ha (around 6563A) of t 6563 ~2.3 
and, therefore, an average extinct i on in the V band of Ay~3mag, using the extinction law 
parametrization by ICardelli et al.l ( 119891 ) and assuming the ratio of total to selective absorp- 
tion, Ry, eq ual to 3.1. Our r esult is in good agre ement with previous estimates of Ay=3.3 
and 4.8 mag (lGoodrichlll99ll ; ICohen &: Kuhilll977l ). especially taking into account the lower 
intrinsic Balmer decrement, 3 and 2.85, adopted in these earlier works. 

We have compared the extiction computed above wi th that predicted from the d etailed 
model of the molecular envelope of M 2-56 performed by ICastro-Carrizo et al.l (120021 ) based 
on high-angular resolution interferometric maps of the CO ( J=l-0 and J=2-l) emission. The 
molecular envelope surrounds the optical nebula and is expected to contain a mixture of gas 
and dust. We have derived the relative variation along the lobes of the CO column density 
from the geometry and density spatial distribution in the model. Adopting a given CO-to- 
H2 relative abundance, X(CO), and the standard conversion factor from H2 column density 
to optical extinction Nh 2 =2.3x10 21 Ay cm -2 , we have calculated the relative variation of 
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the circumstellar extinction (Ay S ) along the lobes produced by the dust in the hourglass- 
shaped molecular envelope (Fig.|9l left). In order to determine the total (circumstellar plus 
interstellar) extinction, we have added the interstellar contribution in the direction to M2-56, 
which is Ay M = 1.75 according to the "Galactic Dust Reddening and Extinction" calculator 
provided by IRSA/IPAC0 (Fig.® right). 



In their model, ICastro-Carrizo et all ( 20021 ) adopt a typical value of X(CO)=2xlO 



-l 



however, these authors affirm that the average CO abundance in this object must be sub- 
stantially lower due to significant photodissociation of this molecule. We find that assuming 
X(CO)~3xlO -5 we are indeed able to reproduce very well the total extinction derived from 
the Balmer decrement. (Note that most of the Ha and, especially, H/3 emission arise in 
the bright ILs, therefore our estimate of Ay~3 mag represents the mean extinction over the 
central ±2"-3"). We have obtained an independent estimate of X(CO) by comparing the 
total dust mass in M2-56 obtained from simple SED modeling, Md us t=l-5xl0~ 3 M Q , with 
the mass of the molecular envelope derived from CO, M gas =O.O5M , assuming the standard 
value of X(CO)=2xl0 -4 . Both masses are related by the following expression: 

M gas = M dust x [^^] x 5, (3) 

where 5 is the gas-to-dust mass ratio, which is typically in the range 5=100-200 in 
pPNs. Applying the previous equation, we deduce an average CO abundance of X(CO)~(3.3- 
6.6) xlO -5 , comparable to that inferred from the extinction and the CO model explained 
above. 



5.1. A compact, dusty "cocoon" around the star? 

We have estimated the total optical extinction along the line-of-sight towards the central 
star of M 2-56 (Ay) from the ratio between the continuum flux observed at a given wavelength 
[F^ hs ) and that expected from the central star in the absence of absorption (Fx)- We 
compute F£ hs from the observed surface brightness at the nebula center and taking into 
account the PSF in our observations. On the other hand, Fx has been approximated by 
the flux (measured on earth) emitted by a black body with an effective temperature of 
T e ff=25,000K, a total bolometric luminosity of L=55OOL , and a distance to M2-56 of 
<i=2.1kpc (§[[]). Considering the values obtained for ^^=1.8 xl0~ 17 ergs -1 cm -2 A -1 and 



http : //irsa . ipac . caltech . edu/applications/DUST/ 
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F65io=4xlO _13 ergs _1 cm -2 A -1 , we infer a value for the total extinction towards the star 
of Ay=12mag. The slope of the continuum emission over the whole wavelength range 
observed by us in 1998 is indeed very well matched by that predicted by the reddened black 
body above. The color of the continuum, however, is not compatible with a lower value of 
Ay ~Ay ~3-4 mag since, in that case, the central star would have to be unreasonably hot 
(T e ff>300,000 K), which is totally in contradiction with the low-excitation conditions and 
lack of substantial photoionization of the optical nebula of M 2-56. 

Our value of Ay is significantly larger than the average nebular extinction implied 
by the Balmer decrement and also larger than the maximum value deduced from the CO 
column density under a reasonable assumption of the CO abundace (see §EJ). This result 
is not totally unexpected: the CO model only accounts for the extinction produced by the 
extended hourglass-like envelope that surrounds the optical lobes, where the Ha and H/3 
emission is locally produced, but does not include the possible attenuation of the stellar 
radiation produced by a potential/conceivable dusty structure inside the optical lobes and 
closer to the stellar surface. The absence of scattered continuum along the lobes in 1998 
suggests that the light from the star was largely attenuated in that epoch not only along the 
line-of-sight, as indicated by the large A v found, but also along the nebular symmetry axis, 
i.e. in the direction to the lobes. It is then quite plausible that the structure responsible for 
the large optical depth towards the star was completely (or almost completely) surrounding 
the latter in 1998 forming, for example, a dusty "cocoon" or shell. Since such a cocoon or 
shell-like obscuring structure is most likely surrounding the close stellar environment and, 
thus, deep inside the ILs, we derive an upper limit to its radius of <0'.'l (<200 AU) from the 
PSF of our high-angular resolution HST images. 



5.2. Recent decrease of the optical depth towards the star 

The remarkable spectral differences observed in 2009 with respect to 1998 and 2000 
(§ !3.2.2p signal important physical changes at the core of M2-56. On the one hand, the 
recent appearance in 2009 of scattered continuum all along the lobes suggests that: 1) the 
intrinsic intensity of the stellar continuum has increased; and/or 2) the optical depth of the 
material (cocoon?) around the star has decreased (at least in the direction to the lobes) 
leading to a larger fraction of the stellar radiation to escape and to reach the lobes where it 
is reflected by the dust. This (option 2) would also explain the increment of the scattered 
(versus locally produced) Ha emission along the lobes in 2009 and, in particular, the recent 
brightening of clumps C3 and C3' evidenced by the direct images (Fig. [I]). On the other 
hand, we have observed a significant brightening of the continuum between 1998 and 2009. 
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In particular, at the nebula center we measure Fggf (2009)~15xFg5f (1998). We believe that 
a true variation of the luminosity or the stellar temperature is very unlikely to be the cause 
of such brightening. First, an increase of the luminosity would be totally unexpected in 
the post-AGB phase, which is characterized by a progressive warming of the central star at 
constant luminosity. Second, only a decrease of the stellar temperature down to T e ff~5500- 
6000 K (at constant luminosity L=55OOL ) combined with an extinction of Ay=11.7mag 
would be able to reproduce the continuum spectrum observed in 2009. Such an extreme 
cooling of the star, however, would not only be unprecedented but, most importantly, would 
be inconsistent with a) the presence of a compact circumstellar H II region in 2009, which 
requires a relatively warm star (T e g >20,000 K), and b) the lack of metallic absorption lines 
typical of a G-type star in the stellar spectrum of M2-56. 

Accordingly, the increase of scattered (continuum and Ha) emission along the lobes and 
the brightening of the continuum at the center are both most likely due to a decrease of the 
optical depth of the circumstellar dust enshrouding the nebula's core (which includes the 
star and the compact circumstellar H II region). 

Applying the same method and adopting the same value for F 6510 as in § 15.11 and taking 
into account the continuum flux measured in 2009 near Ha, F^ =2.7x 10~ 16 erg s -1 cm -2 A -1 
(§ !3.2.2p we deduce a value of Ay=8.8mag. Therefore we conclude that the extinction to- 
wards the star along the line-of-sight has decreased by ~3 mag from 1998 to 2009. 

We have also estimated a lower limit to the decrease of the optical depth from the 
star to the lobes, in particular, to the location of condensation C3 of the eOL (Figs. [I] 
and [5]) from the ratio of the surface brightness of the scattered continuum at that point 
in 2009 and in 1998. The stellar scattered continuum observed at C3 is attenuated by i) 
the dust inside the lobes, i.e. from the star to C3 (symbolized as * — y C3), and also by 
ii) the dust between C3 and us along the line-of-sight, i.e. contained in the extended CO 
envelope that surrounds the optical lobes and in the ISM. However, the extinction produced 
by the latter components is not expected to vary with time (at least, not in less than 
~10yr). In 1998, the continuum is not detected, therefore we use an upper limit given by 
the rms noise in our spectra (lcr=3.1 x 10~ 18 ergs -1 cm -2 A -1 pix -1 ). In 2009, we measure 
1.1 xlO -17 ergs -1 cm -2 A -1 pix -1 at 6510A. Therefore, we conclude that the extinction from 
the star in the direction to the lobes has decreased by Ay* C3 >1.6 mag from 1998 to 2009. 
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6. Density and Nebular Mass 

6.1. The electron density distribution 

We have used the [S II]A6716/A6731 doublet ratio to estimate the electron density, n e , 
along the shocked lobes of M 2-56. In Fig.[10l we plot the spatial profiles along the nebula axis 
(slit PA90) for the two doublet lines and their ratio obtained from the spectrum observed in 
2009 (with the highest S/N ratio and best spatial resolution amongst our multi-epoch data). 
We have checked that the ratios measured in 1998 and 2000 (not shown here) are roughly 
consistent with those obtained in 2009 within the observational errors. 

The [S II]A6716/A6731 ratio increases with the axial distance from the nebula center, 
which implies a decrease of n e from the center to the ou ter parts of the lobes, consid ering 



the [S II]A6716/A6731 ratio versus n e relationship (e.g. lOsterbrock fc Ferlandl 120061 ) and 
adopting an average electron temperature T a = 10,000 K consistent with the estimate by 
Cohen fc Kuhil ( 119771 ) ; iTrammell et al.l (Il993[ ) . We derive the largest value of n e >2 x 10 4 cm -3 



towards the nebula center (within the PSF, i.e. ±0'/5), where the high-density limit of the 
[S II]A6716/A6731 ratio, ~0.4, is reached. Along the ILs, n e ranges between 2000 cm -3 , at 
the base (±0'/5), and 1000 cm -3 , at the tips (±2"). In the inner regions of the OLs (from 
±2" to ±3"), the density steeply varies between <1000 and 500 cm -3 and then decreases 
down to 150 cm -3 at ±5". Beyond these regions the OLs may become more tenuous, if the 
spatially decreasing trend of n e is maintained, however, the large data errorbars prevent an 
accurate determination of the doublet ratio. 

The p-v diagram of the two lines of the [S II] doublet and their intensity ratio, which 
are represented in Fig.[10l provides more details on the distribution of n e in M2-56. In 
particular, unveils the presence of a dense equatorial region at the center (approximately 
within ±1") with a velocity gradient markedly different from (opposite to) that in the lobes. 
Such a velocity gradient is indicative of expansive motions in the direction perpendicular 
to the nebula symmetry axis. The velocity and the density along the equatorial flow are 
both found to increase with the distance to the center, reaching values of up to ~ 100 km s _1 
and >2xl0 4 cm -3 . (We cannot rule out a radial decrease of the electron temperature being 
totally or partially responsible for the radial decrease of the [S II]6716/6731 ratio). The 
center of the expansion may be slightly bluewards of V sys , at VJ.~— 60 km s -1 . 

Finally, we have estimated a lower limit to the density of the region/wind responsi- 
ble for the transient feature Fl observed in the Ha profile (and, to a lesser extent, the 
[N II]AA6548,6584A doublet) in 1998. We believe that the most likely reason why feature Fl 
had vanished in 2000 is that the material ejected during the Ha-emission burst had become 
fully (or almost fully) recombined in the time span between our runs #1 and #2, i.e. in 
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just 2.9 yr or less. This enables estimating a lower limit to the average density in the fast 
Fl-wind in 1998 of n e >2.6xl0 4 cm -3 us ing the expre ssion for the recombination time scale, 
t r , as a function of the density given by Kwokl ( I2000l ): t r ~ 7.6 x 10 4 /[n e / cm -3 ] yr. Such a 
value of the density, relatively high, is consistent with feature Fl not being present in the 
[S II] spectrum: given the critical density of the lines in the doublet these transitions would 
be significantly suppressed by collisional de-excitation in high-density regions. The lines of 
the [N II] doublet have larger critical densities, n c ~fewxl0 4 cm~ 3 , that is why the Fl-wind 
is traced by these transitions. 



6.2. Atomic and Ionized Mass 



We have estimated the total mass of atomic and ionized gas, Mh and Mh+ , in the lobes 
of M2-56 using the mean electron densities derived above and the total energy radiated 
by [O l]A6300 and Ha, respectively. The [O I] (Ha) intensity is proportional to the prod- 
uct of n e and the H (H + ) number density, assuming that the transitions are optically thin 
and that the electron temperature does not strongly vary within the emitti ng region. Con- 
sider i ng a mean electron te mperature in the shocked lobes of T C =10,000K (I Cohen &: Kuhi 
19771 ; iTrammell et al.lll993[ ). relative abundances of He/H = 0.1 and O /H=4xl0~ 4 and the 
required atomic parameters (see e.g. lMendozalll983l ; lGurzadyanlll997l ) we derive: 



Mjj OI] (M ) = 7.7 x 1(T 5 ( 
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(5) 



where L is the dereddened luminosity of the line (and we have adopted d=2100pc). Note 
that Eq. is valid for electron densities smaller than the critical density of [O I], n c ~ 
10 6 cm -3 , which is the case for M2-56. On the other hand, this equation only provides a 
lower limit of the total atomic mass (Mh, which includes the He contribution), since we 
are assuming that most of the oxygen is neutral and in the ground state. Also, we note 
the dependence of the mass derived using Eq. (j3J on T c as M^ 1 ' oc \/{T e ) / '( e ~ 2 - 3e4 / 'T e ) . In 
deriving Eq. ((5]) we have assumed that the majority of the electrons in the nebula come 
from H + , i.e., n e ~ Uh+, and the classic radiative recombination case b for T e ~ 10,000 K 
(i.e. Zx=n ? ./n P nu + =0.25 xlO" 20 cm 3 , where n 3 is the population of the level n = 3 of H; 



Gurzadyanl 1 1 99 71 ) . 
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We have separately estimated the masses of the OLs and ILs using the line fluxes 
integrated within each component and using Eqs. (jlj) and (jSJ). We have also estimated M H + 
in the Fl-wind calculating the flux in the feature Fl computed as the difference between 
the Ha fluxes measured in 1998 and 2000. To derive the total mass in the extended OLs we 
have used the fluxes in Table [3] multiplied by a factor ~3, which converts the flux within a 
l"-wide slit to total nebular flux (estimated from the direct HST and ground-based images, 
FigJT]). We assume, however, that the compact ILs and wind-Fl lie completely within our 1"- 
wide slit. Our results are given in Tab le[3] together with the average values of the extinction 
and electron densities used for each nebular component (estimated in §|5] and § 16. ip . For 
the Fl-wind our value of the ionized mass is particularly uncertain since the density and 
the extinction are poorly known in this case. The extinction by which the flux in the Fl 
feature needs to be corrected for could be of up to ~12, if it is extinguished similarly to 
the central star, or ~3, if it is extinguished similarly to the ILs. Given the angular size 
of the Fl-emitting region (0'.'5), which is larger than that inferred for the cocoon around 
the nebula's core (<0'/l, § 15. ip . we favor a moderate extinction (i.e. closer to ~3) for this 
component. Given the lower limit to n e , we calculate an upper limit to the mass of the 
Fl-wind of <4xlO~ 6 M Q . 

We have also estimated the mass of ionized material in M 2-56 geometrically, using the 
mean electron densities and the volume of the shocked lobes directly measured from the im- 
ages, assuming a simplistic spherical geometry for the lobes. We have separately considered 
the case of hollow lobes and lobes uniformly filled with gas. Assuming hollow lobes with 
~0'.'2-thick walls, the ionized mass of the ILs and OLs is ~7xlO _5 M and ~6xlO~ 4 M , 
respectively. These values are consistent with those previously obtained from the Ha lu- 
minosity (Table[3]). For the adopted values of the mean electron density and taking into 
account the limitations of our necessarily simplistic assumptions made to derive the mass, 
the thin-walled lobes scenario seems to be a better approximation to the structure of the 
lobes of M 2-56 than the filled-lobes geometry: for the latter, the derived masses systemati- 
cally exceed those obtained from the Ha luminosity. We note that the limb-bright nening of 
the ILs and, to a lesser extent, the OLs observed in the images support the hypothesis of 
relatively thin-walled lobes. 



7. Discussion: formation and evolution of M 2-56 

We interpret the nebular morphology and kinematics of M2-56 in the context of the 
AGB-to-post-AGB wind interaction scenario, i.e., resulting from the hydrodynamical in- 
teraction between fast, collimated (jet-like) post-AGB ejections and the slowly expanding 
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envelope expelled during the previous AGB phase (see §[T]). The remarkable "lobe inside 
a lobe" appearance of M2-56 together with the discovery of the recent Ha emission burst, 
and the revelation of other very small-scale structures at the nebula's core points to multiple 
episodes of mass outflow during the post-AGB phase. 

The earliest post-AGB mass-loss event in M2-56 (that we have knowledge of from 
existing data) is probably the one that shaped and accelerated the ex tended hourglass-like 



molec ular outflow enclosing the optical nebula (§HD- As discussed by ICastro-Carrizo et al. 



(120021 ). such an interact ion started approximat ely ~1400yr ago and took place quickly, in 



less than 300 yr (see also lBujarrabal et al.ll200ll ). The smaller kinematical ages of the optical 



OLs and ILs derived from our long-slit spectra, tkin~380 and ~40yr, respectively, suggest 
that the latter resulted from two consecutive and more recent post-AGB ejections. The 
duration of such ejections cannot be asserted based on our data but it may be a small 
fraction of their kinematical ages as in the case of the fast molecular outflow. We find that 
the caps of the optical OLs, which are located inmediately behind of two diametrally oposed 
dense axial clumps traced by CO, show a different kinematics than the main body of the 
OLs. In particular, the OLs-caps are characterized by expansion velocities close to those 
of the aforementioned CO condensations. The kinematical age derived for the OLs-caps, 
ikm~1200yr, is thus similar to that of the CO axial clumps, which may suggest that both 
were shaped and accelerated in the same wind interaction process between the primer post- 
AGB bipolar flow(s) and the AGB envelope. In this scenario, the massive molecular envelope 
would be mainly composed of shocked- AGB wind whereas the OLs-caps could be a remnant 
of the primer post-AGB wind. Alternatively, the material in the OLs-caps could have been 
ejected at the same time as the rest of the OLs, i.e. ~380 yr ago, and subsequently decelerated 
(from expansion speeds of ~400km s _1 down to the observed velocities in these regions, 
V r exp ~110km s" 1 ) as a consequence of its interaction with the pre-existing CO clumps. The 
fast ILs probably trace a more recent shock interaction between an ensuing, fast post-AGB 
wind with the material at the base of the pre-existing and slower OLs. 

Our data also uncover a dense equatorial structure at the nebula center that is expanding 
in the direction perpendicular to the lobes with radial velocities of up to ~100km s~ x (§ 16.11 
and Fig. [TO]) . The velocity gradient observed, which is approximately linear, indicates a 
kinematical age of ~300-400yr (adopting i=90° — 18°=72°). The similar age of the equatorial 
flow and the OLs is consistent with both components resulting from a single, sudden mass- 
ejection event that took place simultaneously along the nebula axis an d in th e perpendicular 



plane. This scenario has been proposed previously by lAlcolea et al.l (120071 ) to explain the 
radial acceleration along the equator and bipolar flows of the molecular envelope of the 
pPN Ml-92. In M2-56, the compact equatorial structure unveiled by the [S II] doublet 
may represent the inner (warmer) regions of the dense, molecular ring/torus probed by CO 
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emission (§[!]). The lower expansion velocity of the latter (\4 xp =7-8km s -1 ) may indicate 
substantial deceleration of the equatorial flow after interaction with the dense, slow AGB 
envelope. The density is found to increase outwards along the optical equatorial flow, in 
support of the presence of compressed (swept up) material at the boundary with the dense 
molecular torus. 

The discovery of the Ha emission burst from the nebular core in 1998 (feature Fl) 
brings to light a latter stellar wind. As discussed in § I4.3[ the Fl-emitting wind is probably 
bipolar. The linear size and mean expansion velocity of this fast, compact flow allow us to 
derive an upper limit to its kinematical age of only <14yr. The upper limit arises because 
of possible deceleration of the Fl-wind by interation with the ILs. The absence of feature 
Fl in the Ha spectrum obtained in 1989 by G91, in fact, suggests that the Fl wind could 
have been ejected only <10yr before it was first identified in our run^l. In this case, the 
initial expansion speed of the Fl-wind (before substantial deceleration took place) should 
have been V^ xp ~500km s _1 to explain the axial extent of this component. This value of V exp 
is indeed in agreement with that inferred from the broad wings and the P-cygni absorption 
component of feature Fl, which indicate that a fraction of the Fl-flow is actually moving at 
this relatively large speeds (§ 14. 3p . 

The absence of feature Fl in 2000 is probably due to fast recombination of the gas 
in the dense (>2.6xl0 4 cm -3 ) Fl-wind (§ 16. ip . In order for the Fl-wind to be completely 
or considerably neutral already in 2000, most of it should have been far away (detached) 
from the central star at that time, otherwise one would expect to observe emission from 
the dense photoionized material at the base of the wind close to the stellar surface as long 
as the wind is on-going. We conclude, therefore, that the rapid recombination observed, in 
less than ~3yr, is most consistent with the Fl-wind being a short-time duration event, e.g. 
a pulsed-jet or bullet-type ejection, leading to detached or clumpy structures, rather than 
being blown in a continuous way. We believe that the clumpy, jet-like structures observed 
in the direct HST images of M2-56 at the base of the elL could represent the footprints 
of such a fast Fl-wind. This is based on the similar size and location of both structures in 
the images and the spectra. If this is correct, the emission that we see from the prominent 
three-finger like structure in the F606W HST image must be mainly stellar light scattered 
by the Fl-wind, after recombination, since that image was obtained in 2002 once the Ha 
emission feature Fl had dissapeared. 

We have estimated a lower limit to the mass-loss rate that led to the main nebular 
components identified in our optical imaging and spectroscopic data, namely, the OLs, the 
ILs and the Fl-wind, from the ratio between the total mass contained in each component 
and their kinematical ages (TableEj). The lower limit arises because the duration of each of 
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the multiple post-AGB mass-loss episodes undergone by M 2-56 could be shorter than their 
kinematical ages. We obtain very similar values of the mass-loss rate for the OLs and ILs, 
namely, Mol s ~2x10~ 6 and Mil s ~3x10 -6 M q yr _1 . For the Fl-wind, we compute a value of 
Mfi~4x1O _7 M yr _1 , however, we note the larger uncertainties in this case mainly due to 
our vague estimate of the mass of the Fl-wind. We have attempted a different estimate of the 
mass-loss rate of the Fl-wind as given by M = 7rr|\4 xp p, where the geometry of the Fl-wind 
is approximated by a cylindrical structure with radius rj, density p, and expanding at V exp . 
Adopting r j ~0'/05, p>3xl0 4 cur 3 and V A exp ~500km s" 1 , we derive M F1 ~4xlO~ 7 M yr _1 . 
The radius adopted for the Fl-wind is that of the bright, jet /finger- like features in the 
F606W HST images that may be footprints of the Fl-wind. The agreement between our 
two estimates of Mpi suggests that the order of magnitude obtained may be correct. 

In general, very little is known about mass-loss rates of post-AGB winds, their con- 
tinuous or episodic nature, and about the structure of the resulting flows after the AGB- 
to- 'post-AGB' wind interaction. Our data indicate that the various post-AGB mass-loss 
episodes experienced by M 2-56 did not happen at regular time intervals, in particular, the 
time span between two consecutive post-AGB ejections has shortened with time. Also, ac- 
cording to the different expansion velocities of the distinct nebular components of M2-56, 
including the molecular envelope, the succesive post-AGB winds seem to be characterized by 
ejection speeds exponentially increasing with time. It is also worth mentioning that, unlike 
most pPNs and PNs with multiple optical lobes, the OLs and the ILs of M 2-56 are oriented 
almost identically (with their main symmetry axis inclined with respect to the line-of-sight 
by i=18°; § !4.ip . Since the molecular envelope has also a similar orientation, we conclude that 
the different post-AGB bipolar ejections in this object (except maybe for the Fl-wind) have 
taken place along the same (or nearly the same) direction/axis over the last ~1400yr. We 
note that, in contra st, signatures of sign i ficant jet directio nal changes are common amongst 



pPNs and PNs (e.g. iMiranda et al.lll999l ; ISahai et al.ll2007l and references therein). A nother 



example of directionally stable jets is found in the PN He 2-90 ( iGuerrero et al.ll200fl ). 



We have obtained a rough estimate of the scalar linear momentum carried by the bipolar 
OLs, the ILs, and the Fl-wind from the product of their masses (Table[3]) by their maximum 
expansion velocities, which have been taken to be ~250, 350, and 500 km s _1 , respectively. 
In this calculation we assume that the nebular components are elongated structures with 
a velocity field that is mainly axial (§ 14.21) . The values found (PoLs~O.17M km s _1 > 
-PiLs~0.05 M km s _1 > Ppi~2xl0~ 3 M km s _1 ) are orders of magnitude smaller than the 
momentum carried by the molecular outflow, P mo i~lOM km s _1 (in deriving this value, the 
linear momentum computed by Castro-Carrizo et al., 2002, P mo i~2 M km s -1 , has been cor- 
rected for our lower estimate of the CO-to-H 2 molecular abundance - see §EJ). According to 
this, we must conclude that none of the post-AGB ejections probed by our optical data could 
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have transferred its large linear momentum to the CO outflow and, therefore, a primer more 
energetic post-AGB ejection is necessary to explain the dynamics of the molecular envelope. 
This result is in good accordance with the time decreasing trend of the linear momentum 
of the optical nebular components inferred by us (Pol s >Pil s >Pfi), which independently 
suggest (by extrapolation) that earlier post-AGB flows could have carried a larger linear 
momentum. 

The presence of a compact dusty structure and a nuclear H II region unveiled by our data 
(§ 15.11 and 14.41) show evidence of further post-AGB mass-loss after the Fl-wind in M2-56. 
The decrease of the optical depth along the line-of-sight and along the nebular symmetry 
axis observed (§ 15.21) is consistent with a detached/shell-like (cocoon?) structure moving 
away from the central star. Assuming that the cocoon expands at constant velocity, v c , and 
that the density varies with the radial distance to the star as p oc r~ 2 (i.e. the mass-loss rate 
is constant), it can be easily demonstrated that for a thick and a thin shell-like geometry: 



where r 98 is the inner and mean radius of the thick and thin cocoon, respectively, in 1998, 
^98 (^09) is the column density in 1998 (2009), and At is the time scale of the variation of 
the optical depth. In 1998 and 2009, the total extinction towards the star was Ay=12mag 
and 8.8 mag, respectively (§ 15. 2p . Subtracting the component of the extinction produced by 
the CO envelope and the ISM at the center, ~3 mag, we estimate that the V^-band extinction 
produced by the cocoon was A^ coon =12-3=9 mag in 1998 and A^ ocoon =8.8-3=5.8 mag in 2009. 
Since the extinction is proportional to the column density, we find iV98/iVo9=9/5.8=1.55. 
Adopting rgg^O'.'l (~3.14xl0 15 cm, § 15. ip and At~10yr, the variation of the optical depth 
observed can be explained for moderate expansion velocities of the cocoon of t> c ~30 and 
60 km s~ x for the thin and thick shell approximation, respectively. 

In the previous expanding cocoon scenario, the detached geometry (needed to explain 
the decrease of the optical depth) implies that the mass-loss process that led to this structure 
had already ended in 2009. The compact size of the cocoon, which is well inside the ILs 
and plausibly closer to the star than the Fl-wind, suggest it is at least as young as the 
latter (i.e., tkm<10yr). Adopting an expansion velocity of t> c =500km s" 1 as for the Fl- 
wind, which is not unreasonable given the increasing trend of the expansion velocity of the 
multiple post-AGB winds of M2-56, we obtain tkm^2yr. 




(thick shell) 



(thin shell), 



(6) 



(7) 
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The presence of compact cocoons enshrouding the central source has been suggested in 
several pPNs. For example, the so-called searchlight beam features that are observed emerg- 
ing from the nebular core in an increasing number of pPNs and PNs (of which CRL 2688 
remains as the best known example) have been h ypothesized to result from starlight escaping 



throu gh holes/cavities in an inner dust cocoon (ISahai et al.lll998l ; ISanchez Contreras et al. 



20071 ). Since these cocoons have not been directly detected, their existence remains specula- 
tive so far. We cannot rule out that the large obscuration of the central source of M 2-56 both 
along the lobes and the equatorial plane could also be produced by two distinct nebular com- 
ponents, namely, a dusty bipolar flow and an equatorial disk/torus, rather than by a unique 
cocoon-like structure. Direct observations of the obscured, innermost regions of M 2-56 us- 
ing, for example, high-angular resolution techniques in the infrarred and mm-wavelength 
range are crucial to characterize the circumstellar geometry at the nebula's heart. 

Finally, our spectra taken in 2009 have revealed a compact H II region around the central 
star of M2-56 that may represent the very latest (current?) mass-loss episode of this object. 
It is possible that such a compact H II region was already formed in 1998 but its light, like 
the stellar radiation, was heavily blocked from our view by the dusty cocoon. Only due to 
the recent decrease of the optical depth of the latter, the emission from the H II region has 
been able to reach the lobes where it is scattered by the nebular dust, enabling its indirect 
detection. Alternatively, it is also possible that the H II region has formed recently (after 
2000). In the following, we discuss this possibility in some detail. The recent formation 
of the H II region would imply an equally recent and rapid evolution (in less than ~10yr) 
of the central star towards higher effective temperatures, needed to explain the raise of 
photonionizing UV radiation and detectable Ha emission from the ionized regions. Although 
a large increase of T e g is not deduced from our multi-epoch spectroscopic data, which do 
not show significant changes in the stellar spectrum^, a low/moderate increase of the stellar 
temperature from a B0.5 to a B0 spectral type (AT off ~3000 K) is simultaneously consistent 
with a stable B-type classificatio n and with a large increase of t he Ha equivalent width from 



W\=50 to 300A (for details, see ISanchez Contreras et al.l 120081 ). The T e g evolutionary rate 



value deduced in this case, AT e g/At>3000/10~300 Kyr 1 , is larger than that expected for 
the central star of M2-56 based on theoretical evolutionary models for post-AG B objects 



taking into account the lo w luminosity and, thus, initial mass of this object (e.g. iBloecker 



19951 ; Ivan Hoof et al.lll997| ). However, these models assume constant values for the post-AGB 
mass- loss rate of 1O~ 7 -1O- 8 M yr 1 , single-star evolution scenarios, and other simplifying 
hypothesis that may not be appropriate for M2-56. In particular, we do not rule out that the 
lattest sudden mass ejections leading to the Fl-wind and the compact cocoon, could have 



6 except for its brightening as a result of the decrease of the cocoon optical depth - see §[3] 
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sped up/boost the evolution of the central star towards slightly higher T e g and the recent 
emergence of a nuclear H II region in just a few years. 

In addition to the central H II region, there are probably shocks currently happening at 
the innermost circumstellar regions of this fascinating object as evidenced by the recent emer- 
gence of Fe II lines from the nucleus (§ !3.2.2p . These Fe II lines are known to be good tracers 
of astrophvsical shocks, e.g ., in supernovae, Herbig Haro objects, PNe, etc (e.g. IWelch et al. 



19991 ; iReipurth et al.l 120001 ) . The presence of shocks at the nebula's core suggests current, 
fast stellar ejections (in agreement with the broad nuclear Ha profile) interacting hydrody- 
namically with the material in the close environment of the star ejected in earlier mass-loss 
episodes with smaller velocities. High-angular resolution optical imaging and spectroscopy 
are needed to properly study the spatio-kinematic structure of the innermost nebular re- 
gions and to attempt understanding the complex mass-loss history of the latest post-AGB 
ejections in M2-56. 



8. Summary 

We report multi-epoch long-slit spectra at various slit positions and direct images in 
the optical of the pPN M2-56. These data have allowed us to describe the spatio-kinematic 
structure and complex (currently active?) mass-loss history of this object with unprecedented 
detail. 



Our data probe several nebular components, namely, two pairs of nested, co-axial 
lobes with different sizes and a large contrast in surface brightness referred to as the 
faint outer lobes (OLs) and the bright inner lobes (ILs), a compact bipolar flow (the 
Fl-wind), and an equatorially expanding central structure. A compact circumstel- 
lar structure obscuring the star (cocoon?) and a nuclear H II region, both spatially 
unresolved, are also inferred from our data. 

The optical lobes are oriented along the East- West direction (PA=90°) and their sym- 
metry axis is inclined ~18° (with respect to the plane of the sky), similarly to the 



hourglass- shaped molecular envelope that surrounds the optical nebula (jCastro-Carrizo et al 



20021 ). 



The OLs and the ILs have an overall axial symmetry and are characterized by an 
expansive kinematics with the velocity increasing with the distance to the nebula cen- 
ter, except for the outermost OLs-caps, which show reduced speeds relative to their 
innermost regions. The maximum velocity reached by the OLs (ILs) is ~250km s _1 
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(~350km s 1 ). Consistent with the large speeds observed at the tips of the ILs, we 
measure proper motions (A^ t ~0'. / 03 yr -1 ) by comparing two-epoch HST images. 

- In 1998, we detected a burst of Ha emission from the nebula nucleus (referred to as 
"feature Fl") that has vanished in less than 2.9 yr. We believe that feature Fl arises 
in a dense, fast (~350-500km s _1 ) bipolar wind ejected after 1989 as a short-time 
mass- loss event. 

- The mass in the OLs, ILs, and Fl- wind is M OLs ~7x 1(T 4 M , Mi Ls ~1.2xl(r 4 M , and 
M F i<4xl(T 6 M , respectively. 

- The p-v distribution of the [S II]6716/ 6731 line ratio has led to the discovery of a dense 
(n e >10 4 cm" 3 ), equatorial flow. The radial velocity increases approximately linearly 
with the distance to the center, reaching values of up to ~100km s -1 . It is possible 
that the equatorial flow and the OLs both resulted from a single, sudden mass-loss 
event that took place simultaneously along the axis and in the perpendicular plane. 

- We infer the presence of a compact (<200 AU), dusty structure (cocoon?) enshrouding 
the central star. The optical depth of this structure has decreased a few magnitudes 
both along the line-of-sight and along the optical lobes in <10yr. This has enabled a 
fraction of the light from the nebula's core to escape the stellar vicinity and reach the 
lobes where it is scattered off by nebular dust, leading to a brightening of the scattered 
stellar continuum and the increase of scattered Ha emission along the lobes in 2009. 

- The scattered Ha emission detected in 2009 along the lobes arises most likely from a 
compact H II region around the star. 

- In 2009, we have discovered a number of permited emission lines arising from the 
nucleus. The detection of Fe II lines suggests the presence of shocks at the stellar 
neighborhood resulting from current stellar wind activity. 

- The kinematical ages of the main nebular components of M 2-56 are different from each 
other: tk~1400yr for the molecular bipolar flow, tk~300-400 yr for the OLs and the 
equatorial flow, t k ~40yr for the ILs, t k <10yr for the Fl-wind, and even smaller for 
the compact cocoon and H II region. 

- The rapid evolution of M2-56 is driven by multipe episodes of mass ejection, that 
is, through a gusty or episodic post-AGB wind, that has led to the nested bipolar 
morphology of the nebula and the younger nuclear components discovered in this work. 
The duration of each of such ejections may be a small fraction of its kinematical age, 
which implies extremely short life-times, of a few years or less, for some of the mass-loss 
episodes in M 2-56. 
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- The various post-AGB mass ejection events experienced by M2-56 did not happen at 
regular time intervals, in particular, the time span between two consecutive post-AGB 
ejections has shortened with time. 

- The successive multiple post-AGB winds in M2-56 are characterized by ejection speeds 
increasing with time from 200 to > 500 km s _1 . In contrast, the mass-loss rate and 
linear momentum may show a time decreasing trend. None of the post-AGB ejections 
probed by our optical data could have transferred its large linear momentum to the 
molecular outflow (P mo iPalOM km s _1 ) and, therefore, a primer more energetic post- 
AGB ejection is necessary to explain the dynamics of the molecular envelope. 

In summary, the data presented here indicate that we are witnessing the current, vari- 
able wind activity at the core of M 2-56 and the rapid evolution of the shocked nebular 
material in very short time scales. This makes of M 2-56 a unique object for studying the 
yet poorly known processes responsible for nebular post-AGB evolution through follow up 
studies. New hydrodynamical symulations of the AGB-to-'post-AGB' wind interaction us- 
ing input parameters for the post-AGB wind in accordance with observational properties 
deduced from this work may be very useful for improving our understanding of PN/pPN 
shaping and evolution. For example, post-AGB winds in the form of a series of sudden, 
short-lived ejection events (not necessarily regularly spaced in time) should be considered. 
Also, models should allow for time variations not only in the velocity of the fast wind but 
also in its mass-loss rate and linear momentum. 
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Fig. 1. — Left) Images of M 2-56 obtained with the WFC of the 2.5 m INT through the nar- 
row Ha and R broad-band filters (top and bottom, respectively) plotted using a logarithmic 
scale. The faint, extended east and west outer lobes (eOL and wOL), nebular condensations 
CI, C2, C3, and C3', and field stars are labeled. Middle) Images obtained with the HST 
through the F656N (Ha) and F606W filters (top and bottom, respectively) plotted using a 
logarithmic scale. The bright, compact east and west inner lobes are labeled (elL and wIL). 
Right) Inset of the HST images (FoV=3'. / 6x3'. / 6) showing the small-scale structure of the 
inner nebular lobes. 
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Fig. 2.— HST images of the ILs of M 2-56 observed in 2002 with the F606W filter (color scale 
and blue contours) and in 1998 with the F656N filter (green contours) showing the proper 
motions of the tips of the lobes and other structural differences within their innermost regions 
(closest to the center). 
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Fig. 3. — Long-slit Ha spectra in M2-56 for the five slit positions observed in our f998 and 
2000 runs. The ground-based Ha image of the nebula is shown on the leftmost panel; slit 
positions are superimposed on the image. The origin of the spatial scale in the spectra coin- 
cides with the point of maximum continuum and Ha emission. The LSR systemic velocity 
of the source, V^^— 27 km s -1 , is indicated by a vertical dashed line on each spectrum. The 
intense, blue-shifted emission feature at the nebula center observed in 1998 and its broad 
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Fig. 4. — Long-slit spectra over the whole spectral range observed in 1998 and 2000 for the 
central slit PA90 (see Tableland Fig.E]). 
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Fig. 5. — Top panel) Long-slit spectrum for slit PA90 obtained in 2009 in a logarithmic color 
scale. The spectrum has been smoothed with a flat-topped rectangular kernel of dimension 
3x3 pixels. The units of the wedge at the top of the box are ergs -1 cm~ 2 A -1 arcsec" 1 . 
Middle and bottom panels) Same as above but continuum substracted and showing 
selected wavelength ranges. Lines detected are labeled (see Table[2]). The blue vertical 
lines are placed at V sys for each of the observed transitions. Contour levels are: 2a, 4a, 
Qa, 8a, 10a, 15a, 20a, 25a, 30a, 40a, 50a, QOa, from 75a to 300cr in steps of 25a, from 
300a to lOOOcr in steps of lOOcr (black), and — 3a,— 6a,— 9a... (magenta) with <7=1.2xlO~ 19 
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Fig. 6. — Left) Superposition of the Ha long-slit spectrum along PA90 obtained in 2009 
(grey color scale) and in 1998 (green contours). The systemic velocity is indicated by the 
vertical dashed line. The loci of clumps C3 and C3' along the slit is shown. Right) One 
dimensional Ha profile extracted within a l"xl" region centered at condensations C3 (blue) 
and C3' (red) and at the nebula center (black). 
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Fig. 7. — Top) Onedimensional spectra from a central region of l"xl" extracted from our 
long-slit spectra PA90 obtained in 1998 (black), 2000 (red), and 2009 (blue). The 1998 and 
2000 spectra have been boxcar smoothed to match a common value of the dispersion of 
~0.7A. Bottom) The same as in top panel after fitting and subtracting the continuum and 
for selected transitions detected in all three epochs using a velocity scale in the x-axis. The 
Ha profile in 2009 has been scaled down by a factor 3.5 for an easier comparison with the 
weaker emission profile in 1998 and 2000. The systemic velocity is indicated by the vertical 
blue-dotted line. Note the similar profile of the [S II] doublet lines in all epochs and the 
remarkable differences in the Ha and [N II] lines. 
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Fig. 8. — Left) Schematic geometry of M2-56 deduced from the model fitting (§ 14. 2p . This 
plot represents a cut of the nebula by a plane perpendicular to the plane of the sky in the 
direction of PA90. In this diagram 1" corresponds to 3.14xl0 16 cm (d=2.1kpc, §[[]). The 
velocity field in the ILs (OLs) is indicated by blue (red) arrows, respectively. Right) Ha 
spectrum along PA90 slit orientation as observed in 1998 and 2000. The position-velocity 
diagrams resulting from the model fitting are superimposed. Contours and color scale are 
the same as in Fig.[3j 




Fig. 9. — Spatial distribution of the circumstellar (left) and total, i.e. circumstel- 
lar+interste llar, extinction (right) alon g PA90 inferred from the model of the molecular 
envelope by ICastro-Carrizo et al.l (120021 ) . The value for the CO-to-H2 abundance ratio, ISM 
extinction, and the conversion factor between H2 column density and dust extinction adopted 
are indicated. 
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Fig. 10.— Top) Spatial distribution along PA90 of the lines of the [S II]AA6716,6732A dou- 
blet (left) and their relative intensity (right) observed in 2009. Bottom) Position-velocity 
diagram along PA90 of the [S II] lines and their ratio observed in 2009. The p-v distribution 
of the 6716/6731 intensity ratio towards the center evidences a dense, equatorially expanding 
structure. 
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Table 1: Log of spectroscopic observations 



Run 


Date 


Telescope 


Grating 


Dispersion 


Range 


Pixel 


Slit 


Exposure 


# 


(yyyy-mm-dd) 


+Instrument 




(Apix- 1 ) 


(A) 


(") 


Position 51 


(s) 


1 


1998-01-17 


INT+IDS 


R1200Y 


0.39 


6370-6770 


0.33 


PA90 


3600 


1 


1998-01-18 


INT+IDS 


R1200Y 


0.39 


6370-6770 


0.33 


PA0,PA0E 


3600,3600 


2 


2000-ll-ll b 


INT+IDS 


R1200Y 


0.22 


6234-6805 


0.19 


PA90,PA0E 


10800,9000 


2 


2000-11-12 


INT+IDS 


R1200Y 


0.22 


6234-6805 


0.19 


PA0W,PA119 


11700,10800 


2 


2000-11-13 


INT+IDS 


R900V 


0.30 


4558-5330 


0.19 


PA90 


27300 


3 


2009-08-23 


NOT+ALFOSC 


VPH#17 


0.26 


6330-6850 


0.19 


PA90 


6800 



a See Fig. [3] for slit position labels. 
& Non-photometric conditions 
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Table 2. Multi-epoch line fluxes observed (F\) and derredened for M2-56 



'10 15 ergs 1 cm 



'10 15 ergs 1 cm 2 ) 



Line 



1998 



2000 



2009 



1998 



2000 



2009 



rip A^toOl.O 










ouyy ) 




[W 11IJ /VOUUU.O 




n n7iTi nzn 






Z.Ul ±.Z J 




[FeII]AA5158.0,5158.8 




0.20(0.04) 






4.9(1.4) 




[N I]A5200.3 




0.65(0.04) 






15(3) 




[Fe II]A5261.6 




0.05(0.04) 






1.1(0.9) 




[O I]A6300.3 




6.4(0.4) 






76(12) 




[O I]A6363.8 




2.0(0.4) 


2.2(0.14) 




23(6) 


25(4) 


Fe IIAAGSSS.S^SSS.S 1 " 






1.12(0.07) 






13(2) 


Fe IIA6416.9 1 " 






0.21(0.07) 






2.4(0.9) 


Fe IIA6432.7 1 " 






0.42(0.07) 






4.7(1.1) 


Fe IlA6456.4t 






0.77(0.07) 






8.5(0.14) 


Co I?A6490.3 t 






0.21(0.07) 






2.3(0.8) 


Fe IlA6516.lt 






0.91(0.07) 






10(2) 


[N II]A6548.0 


2.3(0.2) 


1.6(0.4) 


4.9(1.1) 


24(4) 


17(5) 


52(14) 


Ha A6562.8 


23(0.4) 


15.1(0.4) 


161(1.1) 


240(40) 


160(6) 


1700(300) 


[N II]A6583.5 


6.7(0.3) 


4.8(0.4) 


10.5(1.1) 


70(11) 


50(2) 


110(20) 


Fe lA6663.4t 






0.14(0.07) 






1.4(0.2) 


[S II]A6716.4 


2.6(0.2) 


2.4(0.4) 


2.6(0.7) 


25(4) 


23(5) 


25(14) 


[S II]A6730.8 


2.9(0.2) 


3.2(0.4) 


3.4(0.7) 


28(5) 


31(6) 


33(14) 



Lines discovered in this work. 



Note. — Errors are given in parenthesis. 
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Table 3: Atomic and ionized mass of M2-56 



Component 


(10 -15 ergs -1 cm -2 ) 


F [0 i] 

(10~ 15 ergs -1 cm -2 ) 


<A V > 
(mag) 


<n e > 
(cm -3 ) 


(M ) 


M [Oij 

(Mq) 


OLs 


5.7(±0.2) 


2.0(±0.6) 


2 


200 


5.3xl0" 4 


1.5xl0 -4 


ILs 


9.3(±0.2) 


5.4(±0.4) 


3 


1500 


8.2xl0 -5 


3.8xl0 -5 


Fl-wind 


8.0(±0.2) 




^3 (±12) 


>2.6xl0 4 


4.1xl0 -6t (<4.2xl0 -3 ) 





' Particularly uncertain. See text in § 16.21 



